JOURNAL OF GEOPHYSICAL RESEARCH: SPACE PHYSICS, VOL. 118, 7654-7664, doi:10.1002/2013JA019372, 2013 


STEREO and Wind observations of intense cyclotron harmonic waves 
at the Earth’s bow shock and inside the magnetosheath 

A. W. Breneman, 1 C. A. Cattell, 1 K. Kersten, 1 A. Paradise, 1 S. Schreiner, 1 P. J. Kellogg, 1 
K. Goetz, 1 and L. B. Wilson III 2 

Received 27 August 2013; revised 12 November 2013; accepted 24 November 2013; published 12 December 2013. 

[i] We present the first observations of electron cyclotron harmonic waves at the 
Earth’s bow shock from STEREO and Wind burst waveform captures. These waves are 
observed at magnetic field gradients at a variety of shock geometries ranging from 
quasi-parallel to nearly perpendicular along with whistler mode waves, ion acoustic 
waves, and electrostatic solitary waves. Large amplitude cyclotron harmonic waveforms 
are also observed in the magnetosheath in association with magnetic field gradients 
convected past the bow shock. Amplitudes of the cyclotron harmonic waves range from a 
few tens to more than 500 mV/m peak-peak. A comparison between the short (15 m) and 
long (100 m) Wind spin plane antennas shows a similar response at low harmonics and a 
stronger response on the short antenna at higher harmonics. This indicates that 
wavelengths are not significantly larger than 1 00 m, consistent with the electron cyclotron 
radius. Waveforms are broadband and polarizations are distinctively comma-shaped with 
significant power both perpendicular and parallel to the magnetic field. Harmonics tend to 
be more prominent in the perpendicular directions. These observations indicate that the 
waves consist of a combination of perpendicular Bernstein waves and field-aligned waves 
without harmonics. A likely source is the electron cyclotron drift instability which is a 
coupling between Bernstein and ion acoustic waves. These waves are the most common 
type of high-frequency wave seen by STEREO during bow shock crossings and 
magnetosheath traversals and our observations suggest that they are an important 
component of the high-frequency turbulent spectrum in these regions. 
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1. Introduction 

[ 2 ] High-frequency plasma waves have been known for 
some time to play an important role in the dissipation of 
bulk flow energy across the bow shock [e.g., Papadopoulos, 
1985]. Initial studies of turbulence at the bow shock [e.g., 
Fredricks et al., 1970] used time-averaged spectral data and 
showed the presence of two strongly enhanced spectral com- 
ponents: electromagnetic waves at <200 Hz and electrostatic 
waves from 200 to 800 Hz. Part of the electromagnetic 
spectrum was identified with fluxgate magnetometer data 
on IMP 6 to be composed of whistler mode waves by 
Fairfield [1974]. The electrostatic component was identi- 
fied by Rodriguez and Gurnett [1975] as ion acoustic waves 
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because of a dominant polarization along the magnetic field. 
This identification was supported by Fuselier and Gurnett 
[1984] who determined that the ion acoustic waves had 
wavelengths on the order of A > ItcXq, where A/j is the 
Debye length, using wave/antenna interference effects that 
occur when the wavelength is of the same order or shorter 
than the antenna length. 

[ 3 ] Waveform measurements have the advantage over 
low cadence spectral data of not averaging over a certain 
time period. This is particularly important when sampling 
regions with strong gradients and hence short-scale lengths, 
like the bow shock. The first results using electric field 
waveforms (ISEE 1) at the bow shock were published by 
Formisano and Torbert [1982] who identified ion acoustic 
waves downstream of the transition region. Later Wygant 
et al. [1987] identified large amplitude electric field spikes 
from 32 Sample/s survey electric field data in the transi- 
tion region, interpreted as Doppler-shifted lower hybrid or 
ion acoustic waves. Initial waveform results from the Wind 
satellite mission showed the existence of large amplitude 
electrostatic solitary waves at the shock transition region 
[Bale et al, 1998]. These structures, which tend to appear as 
small amplitude broadband bursts in time-averaged spectral 
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Table 1 . Data on All Burst Waveform Captures (# Captures) Observed During STEREO Spacecraft (SC) Bow Shock Crossings 3 


Date SC 

Time 

V 
’ sw 

# Captures 

Sample rate 

Length 

$Bn 

|Emax | 

l^minl 

|Emean | 

|Emedian | 

GSE Position 

06 Nov A 

00 : 51 

380 

17(13) 

125 

32.7 

90 

273(273) 

6(37) 

50(88) 

38(62) 

[15.1,1.8,5.8] 

06 Nov B 

00 : 18 

380 

22(19) 

125 

32.7 

79 

357(357) 

31(163) 

194(268) 

195(265) 

[13.9,1.9,5.6] 


14 : 48 

370 

24(17) 

125 

32.7 

30 

364(364) 

16(109) 

118(176) 

109(163) 

[11. 8-7.2, -4.7] 

17 Nov A 

05 : 51 

390 

38(19) 

250 

16.4 

58 

381(381) 

9(50) 

95(197) 

79(191) 

[14.1-0.9,6.5] 

18 Nov A 

00 : 40 

340 

14(4) 

250 

16.4 

85 

157(137) 

2(67) 

41(95) 

26(95) 

[14.4-11.0,-6.6] 

17 Nov B 

02 : 44 

395 

5(0) 

250 

16.4 

47 

121(N/A) 

6(N/A) 

32(N/A) 

27(N/A) 

[15.2-0.9,6.3] 


20 : 13 

340 

38(22) 

250 

16.4 

20 

225(225) 

4(57) 

58(116) 

49(108) 

|12.2,-10.3,-5.7| 

29 Nov A 

11 : 16 

460 

10(0) 

125 

131 

46 

263(N/A) 

4(N/A) 

52(N/A) 

36(N/A) 

[14.6-2.0,7.0] 

30 Nov A 

03 : 58 

400 

2(2) 

31 

2100 

27 

518(518) 

414(414) 

466(466) 

466(466) 

[8. 1,-0. 1,5. 3] 


04 : 03 

400 

3(3) 

31 

2100 

27 

444(444) 

126(126) 

324(324) 

400(400) 

[8. 1,-0. 1,5. 3] 


04 : 15 

405 

3(3) 

31 

2100 

27 

581(581) 

371(371) 

447(447) 

424(424) 

[8. 1,-0. 1,5.3] 

29 Nov B 

10 : 59 

450 

6(0) 

125 

131 

24 

249(N/A) 

21(N/A) 

83(N/A) 

73(N/A) 

[14.5, -1.5, 6.4] 

30 Nov B 

03 : 55 

460 

2(1) 

31 

2100 

27 

520(520) 

325(520) 

422(520) 

422(520) 

[8. 1,-0. 1,5. 3] 


a |£ ma x| an d | Emin | represent the maximum and minimum wave amplitudes (mV/m peak-peak) out of the set of burst captures at each shock. The mean 
|E me an| and median |£ me dian| values are calculated from all of the burst captures at each shock. The values in the parentheses are for the subset of waves 
with clear comma-shaped polarization only. Solar wind velocity (T sw ) are from the ACE spacecraft. Shock normal angles (#Bn) are estimated from the 
geometric model described in the text. Also listed are the burst sample rate (kSamples/s), length of each individual burst (ms), and GSE position (in units 
of Earth radii) of the spacecraft during the time of the observations. 


data, modify the shock potential structure and can provide 
significant energy dissipation at the shock transition region. 
Waveform data from Polar \ I lull et al., 2006] and Cluster 
[Giagkiozis etal, 2011] have also been used to identify large 
amplitude (~80 mV/m) ion acoustic waves at the bow shock 
transition region. An important observation from these and 
other waveform analyses [Wilson et al., 2007, 2010] is that 
these waves can be large amplitude and bursty. The com- 
bined influence of significant numbers of large amplitude 
waves can potentially provide an important contribution to 
overall energy dissipation at the bow shock. 

[ 4 ] We present a survey of STEREO and Wind burst 
captures showing that, in addition to whistler mode radia- 
tion, ion acoustic waves, and electrostatic solitary waves, 
a significant portion of the electric field spectrum at the 
bow shock transition region consists of electron cyclotron 
harmonic waves. These waves are also observed down- 
stream of the bow shock in the magnetosheath in association 
with strong magnetic field gradients. These waves are much 
larger in amplitude than any previously reported in the 
magnetosheath. 

2. Instrumentation and Data Set 

[5] STEREO observations come from 13 inbound and 
outbound passes through the bow shock as well as a sin- 
gle magnetosheath traversal, all during the phasing orbits 
early in the mission. STEREO Waves [Bougeret et al., 2008] 
measures electric field from three 6 m cylindrical antennas 
[Bale et al, 2008]. These antennas were designed to pro- 
vide accurate measurement of the electric field at frequencies 
> 20 kHz. The accuracy of the electric field measurements at 
the frequencies of waves presented in this paper (~1 kHz) is 
not as well established but it is likely that the amplitudes we 
present are close to their actual values (see Appendix A for 
details of the process of converting the raw antenna signal 
into an electric field). 

[6] The Time Domain Sampler (TDS) provides short 
bursts of AC-coupled electric field at up to 250,000 samples 
per second. No search-coil magnetic field data are avail- 
able on STEREO. The primary advantage of the STEREO 
data set is in providing a large number of three-dimensional 
burst captures during the shock crossings. Shock crossings 


are identified with a combination of In-situ Measurements 
of Particles and Coronal mass ejection Transients (IMPACT) 
fluxgate magnetometer data [Luhmann et al., 2008] and 
Block Adaptive-tree Solar-wind Roe-type Upwind Scheme 
(BATS-R-US) MHD simulations (with solar wind input) of 
the interaction between the solar wind and the Earth’s mag- 
netosphere. In general, particle data are not available to 
facilitate this identification because the STEREO spacecraft 
were still in the commissioning phase during this time and 
the instruments were not fully operational. Despite this, flux 
levels are high enough on a few of the crossings that data 
from the Solar Wind Electron Analyzer (SWEA) | Sam and 
et al, 2008] were used to corroborate the identification of 
the bow shock crossing. 

[7] Wind observations come from a single bow 
shock crossing, identified from fluxgate magnetometer 
[Lepping et al., 1995] and plasma instrumentation [Ogilvie 
et al, 1995]. Wind measures electric field from a short spin- 
axis 12 m tip-tip antenna pair ( E z ) and spin plane antenna 
pairs with tip-tip lengths of 100 m (E x ) and 15 m (E v ). Note 
that on 3 August 2000 and 25 September 2002, the long 
spin plane antenna was sheared by dust impacts, reducing 
its physical length (this effect is accounted for, as discussed 
in detail toward the end of Appendix A). Wind fast mode 
Time Domain Sampler measures short AC-coupled bursts 
of the spin plane electric field [Bougeret et al, 1995] at 
125,000 Samples/s. The different response on the long and 
short Wind spin plane antennas as a function of harmonic 
number allows us to place limits on the wavelength of the 
cyclotron harmonic waves. 

| x | We also show evidence in Appendix B that the har- 
monics observed are real and not introduced by a harmonic 
instrumental response due to large amplitude or short wave- 
length waves. 

3. STEREO and Wind Survey and Event Analysis 

[ 9 ] Table 1 records shock properties as well as burst cap- 
ture statistics at 13 STEREO bow shock crossings. Because 
no particle measurements are available, a rough estimate of 
the shock normal angle at the location of each crossing was 
determined by inputting magnetometer data into the hyper- 
bolic shock model of Slavin and Holzer [1981]. All of the 
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Table 2. Data on All Burst Waveform Captures (# Captures) Observed by the STEREO Spacecraft (SC) Inside of the Magnetosheath 3 


Date SC 

Start Time 

End Time 

# Captures 

Sample Rate 

Length 

l^maxl 

[Eminl 

l-^mean | 

l^median | 

GSE Position Start 

GSE Position End 

17 Nov A 

21 : 54 

23 : 35 

21(21) 

250 

16.4 

110(110) 

45(45) 

67(67) 

65(65) 

[11.3-9.3, -6.1] 

[13.6, 10.6-6.5] 

17 Nov B 

04 : 07 

04 : 38 

21(6) 

250 

16.4 

381(141) 

6(16) 

169(52) 

153(35) 

[13.3-0.4,6.0] 

[13.4-0.4,6.0] 

30 Nov A 

00 : 05 

00 : 42 

9(7) 

31 

520 

280(260) 

110(110) 

183(179) 

180(180) 

[3. 2, -7. 7,-5. 0] 

[4. 0,-8. 4,-5. 2] 

30 Nov B 

00 : 03 

00 : 04 

2(2) 

31 

520 

180(180) 

90(90) 

135(135) 

135(135) 

[3. 8, -8. 8,-4. 7] 

[3.8,-8.8,-4.7] 


“l^maxl and 7/ mm represent the maximum and minimum wave amplitudes (mV/m peak-peak) out of the set of burst captures in each row. The mean |is me an 
and median |£ lnel j iall | values are calculated from all of the burst captures in each row. The values in the parentheses are for the subset of waves with clear 
comma-shaped polarization only. Also listed are the burst sample rate (kSamples/s), length of each individual burst (ms), and GSE start and end positions 
(in units of Earth radii) of the spacecraft during the time of the observations. 


shock crossings concurrent with burst waveform captures 
occur during times of slow to moderate solar wind speed. 
Shock normal angles vary from quasi-parallel (STEREO- 
B on 6 November 2006, 14:48) to nearly perpendicular 
(STEREO-A on 6 November 2006 at 00:51). Of note is that 
large amplitude electric field burst waveforms, often up to 
a few hundred mV/m, are seen at all 13 listed shock cross- 
ings, regardless of shock normal angle. A majority of these 
bursts occur at the shock ramp (when a clear ramp is identifi- 
able). Though the waves are bursty, their amplitudes usually 
far exceed the typical few tens of mV/m quasi-static elec- 
tric field component [Eastwood et al., 2007; Dimmock et al., 
2012 ], 

[ 10 ] Four types of waveforms are seen during the bow 
shock crossings: electrostatic solitary waves, whistler mode 
waves, ion acoustic waves, and cyclotron harmonic waves. 
A small number of the burst captures show clear evidence 
of electrostatic solitary waves as reported by Bale et al. 
[1998]. These structures are usually immersed in a wavefield 
of a few mV/m. Well-isolated solitary waves are not com- 
monly observed in the STEREO burst data set. Nearly all of 
the burst captures show at least some power at frequencies 
/ < fee, where f ce is the electron cyclotron frequency, which, 
when elliptically polarized, suggests whistler mode waves. 
Occasionally, the whistler mode signal dominates the burst 
capture. However, a majority of the burst captures consist 
of large amplitude nonlinear fluctuations lasting for tens of 
milliseconds (few tens of wave cycles) with power peaked 
at / > f ce . It is likely that some of this power is attributable 
to short wavelength ion acoustic waves Doppler-shifted in 
the spacecraft frame to higher frequencies. These show no 
tendency to be polarized along the magnetic field, consistent 
with observations by Fuselier and Gurnett [1984], Balikhin 
et al. [2005], and Hull et al. [2006]. 

[n] A majority of burst captures with peak power at 
f > fee cannot be attributed solely to ion acoustic waves. 
These waveforms have broadband power spectra above f ce 
with significant power both parallel and perpendicular to 
the magnetic field resulting in distinctive “comma-shaped” 
polarizations. These waveforms are the focus of our study 
and their statistics are compiled separately in Table 1. A 
subset of these burst captures shows unambiguous evidence 
of multiple cyclotron harmonics and are identified as elec- 
tron cyclotron harmonic waves. These are the first reported 
observations of electron cyclotron harmonic waves at the 
bow shock. 

[ 12 ] Not listed in Table 1 are the four bow shock cross- 
ings on 12 December 2006 (two on each spacecraft) where 
no burst waveform captures were telemetered. STEREO 
TDSMax, a once-per-minute record of the maximum 


high-frequency electric field, indicates that large amplitude 
waves were indeed present on all four crossings. The lack 
of burst captures then suggests that they were overwrit- 
ten by later large amplitude waveforms, like the extremely 
large amplitude whistler mode waves reported by Cattell 
et al. [2008], before they had a chance to be telemetered 
to the ground. 

[ 13 ] Large amplitude waveforms with comma-shaped 
polarizations, enhanced broadband power at f ce , and 
cyclotron harmonics were also observed on STEREO in 
the magnetosheath. Table 2 presents statistics of these 
burst captures. 

[ 14 ] In the next sections we present a detailed analy- 
sis of burst captures with distinct electron cyclotron har- 
monics at two different bow shock crossings and a burst 
capture associated with a magnetic ramp-like structure in 
the magnetosheath, likely convected downstream from the 
bow shock. 

3.1. STEREO Bow Shock Example 

[ 15 ] We now discuss electron cyclotron harmonic waves 
observed at a bow shock crossing on STEREO-B and shown 
in Figure 1. Figure la plots IMPACT magnetometer data 
showing a bow shock crossing at 20:14 and a bow shock 
crossing or upstream structure from 20:28 to 20:40. Also 
plotted are the locations and amplitudes of the burst captures. 
All of the burst captures are clearly associated with mag- 
netic field gradients. Each shows significant wave power at 
/ > fee, with amplitudes up to 300 mV/m peak-peak. A few 
of the waves at and near the large magnetic field gradient 
at 20:14 show unambiguous cyclotron harmonics. Figure lb 
shows a typical example of a waveform with cyclotron har- 
monics in a magnetic field-aligned coordinate system where 
the magnetic field points along z and the maximum variance 
direction lies in the x-z plane. The waveforms are irregu- 
lar and asymmetrical giving them a “choppy” quality. The 
spectra in Figure lc show that the harmonics exist primar- 
ily perpendicular to the magnetic field. The parallel direction 
also contains significant wave power but without clear har- 
monics. Wave hodograms are presented in Figure Id and 
have the distinctive comma shape in the maximum variance 
plane x-z. As indicated in Table 1, 21 of 38 burst captures at 
this bow shock crossing have comma-shaped polarizations. 
The polarizations for the remaining waveforms are linear 
(without an overall preference for orientation to the magnetic 
field, shock normal, or solar wind velocity) and complex (no 
preferred direction). 

[ 16 ] For the example shown in Figure Id, the comma- 
shaped polarization combined with the observation that the 
harmonics are perpendicular to the magnetic field suggest an 
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Figure 1. (a) STEREO-B observations at a bow shock crossing. The diamonds show the locations and 
amplitudes of the burst captures, many of which show clear signs of cyclotron harmonics, (b) Example 
burst waveform capture indicated by the red vertical line in Figure la. The waveform is in field-aligned 
coordinates where z is the magnetic field direction and the maximum variance electric field is in the x-z 
plane, (c) Spectra of the burst waveform components showing clear cyclotron harmonics (black vertical 
lines) in the directions perpendicular to the magnetic field (black and green) but an absence of harmonics 
along the magnetic field, (d) Waveform hodograms indicating significant wave power in directions both 
parallel and perpendicular to the magnetic field. 


interaction between two types of waves, a Bernstein wave 
perpendicular to the magnetic field and a field-aligned wave, 
possibly ion acoustic. Another possible explanation for the 
comma-shaped polarizations is waveform distortion due to 
electron trapping in large amplitude electric field waves as 
described by Kellogg et al. [2010, 2011]. 

3.2. Wind Bow Shock Example 

[17] Our next example comes from Wind burst data dur- 
ing an inbound quasi-perpendicular bow shock crossing on 
14 November 2001. Similar to Figure la, Figure 2a shows 
fluxgate magnetometer data of the shock crossing as well as 
the locations and amplitudes of the 11 burst captures, 7 at 
the magnetic ramp, with amplitudes (spin plane components 
only) from 80 to 510 mV/m peak-peak. Out of the these 11 
burst captures, one contains discrete solitary waves and the 
rest have clear enhancements of power at or above f ce with 
evidence of cyclotron harmonics. 

[is] Figure 2b shows an 80 mV/m peak-peak nonlinear 
waveform in the Wind spin plane at the ramp. The wavelet 
transform shows that the power is strongly peaked at f ce (as 
determined from Wind high-resolution magnetometer data). 
The spectra in Figure 2c show a clear pattern of harmon- 
ics on both antennas at nearly exact multiples of the local 


value of f ce , particularly on the E y antenna. The magnetic 
field at this time is 7 /gse = [10,-10,9] nT and is inclined 
from the spin axis E x antenna by ~42° and the E y antenna 
by ~90°. Thus, the harmonics are most distinct perpendicu- 
lar to the magnetic field, similar to the STEREO example in 
Figure 1. As seen from Figure 2a, the wave is observed in 
a turbulent region where the magnetic field strength rapidly 
changes over a short distance. Thus, the close association 
of the first harmonic and f ce at the time of observation sug- 
gests that the wave is being observed close to its source. 
In addition, because the visible harmonics strongly peak at 
multiples of //,,, there is likely little Doppler-shift, as would 
be observed if the wave propagated perpendicular to the 
solar wind flow. Three-dimensional waveform data are not 
available to explicitly show this but Figure 2d shows that the 
spin plane polarization is indeed primarily perpendicular to 
the solar wind velocity. 

[ 19 ] The spectra in Figure 2c indicate that the response 
on the two Wind spin plane antennas is nearly identical at 
f ce but that the short (E v ) antenna has a stronger response at 
higher harmonics. This suggests two things: that the wave- 
length of the component f = f ce is longer than the long 
( E x ) antenna length of 100 m or else the response for this 
frequency would be stronger on the short E r antenna, and 
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Figure 2. An example of a cyclotron harmonic wave in the context of a Wind bow shock crossing on 
14 November 2001. (a) Wind magnetic fields investigation (MFI) high-resolution definitive data with 
overlain lower cadence magnetic field data in red to more clearly show the ramp region. The locations 
of all the burst captures are shown with the diamonds with peak-peak amplitudes given by the scale to 
the right, (b) The spin-plane electric field waveform components of the cyclotron harmonic wave at the 
magnetic ramp, indicated in Figure 2a by the vertical red line, along with the wavelet transform of the 
E r component, (c) Fast Fourier transform (FFT) of the waveform showing a peak in power at the local 
value of f ce as well as the first 10 harmonics (vertical black lines), (d) The spin-plane polarization of 
the waveform in Figure 2b. Wind antenna coordinates are shown by the axes labeled E x and E y . Also 
shown are the solar wind velocity and magnetic field vectors. The spin-plane component is predominantly 
perpendicular to the solar wind velocity. 


that the higher harmonics have wavelengths that are on the 
order of or shorter than 100 m. This estimated wavelength 
is roughly consistent with the electron cyclotron radius 
of ~500 m. 

3.3. STEREO Magnetosheath Example 

[ 20 ] We conclude this section with an analysis of a 
cyclotron harmonic wave observed by STEREO-B down- 
stream of the bow shock in the magnetosheath at a mag- 
netic field structure. Figure 3a shows that this structure lies 
between the bow shock and magnetopause. As with other 
large-scale magnetic structures in this region, it likely has 
its origin at the bow shock or upstream of it. Figure 3b 
is a close-up of this structure and shows the locations and 
amplitudes of the burst captures which are up to 100 mV/m 
peak-peak. Like those observed at the bow shock cross- 
ings, these waves are strongly associated with magnetic field 
enhancements and gradients. A majority of these burst wave- 
forms have the distinctive comma-shaped polarization and 
show clear cyclotron harmonics, often most prominently 
in the directions perpendicular to the magnetic field. Past 


studies have shown electron cyclotron harmonic waves to be 
common both inside of the magnetosphere and near the mag- 
netopause but fewer observations have been made within the 
magnetosheath near the bow shock [e.g., Anderson et al., 
1982] and these waves were not reported to be nearly as 
large in amplitude or to have any correlation with magnetic 
field gradients. 

[ 21 ] Figures 3c-3e show an example of one of these 
waveforms in the same field-aligned coordinate system 
described previously. The irregular, asymmetric waveforms 
(Figure 3c) are similar to those seen at the bow shock which 
show signs of electron cyclotron harmonics. The FFTs in 
Figure 3d as well as the hodograms in Figure 3e are also very 
similar to those from waves at the bow shock. Namely, the 
power peaks just above f ce , evidence of harmonics is most 
distinct in directions perpendicular to the magnetic field, 
and the hodograms show the typical comma-shaped polar- 
ization. Again, these results indicate that these waveforms 
are likely a combination of a perpendicular Bernstein wave 
and a parallel wave without harmonics. In this particular 
example there is significantly less power in the field-aligned 
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Figure 3. Example of a cyclotron harmonic wave observed on STEREO-B at a magnetic field structure 
inside of the magnetosheath, (a) The location of the structure relative to the bow shock and magnetopause 
is shown, (b) A closer look at the magnetic structure in the magnetosheath. The diamonds show the 
locations and amplitudes of the burst captures, (c) The burst waveform corresponding to the vertical red 
line in Figure 3b in field-aligned coordinates where z is the magnetic field direction and the maximum 
variance electric field is in the x-z plane, (d) FFTs of the waveform. The local electron cyclotron frequency 
is indicated with the vertical black line. Evidence of harmonics can be seen in all three components but 
are most evident in the perpendicular directions, (e) Hodograms showing the comma-shape that is often 
associated with these waveforms. 


direction, though this may be due to a local field rotation in 
this region of strong field gradient. 

[ 22 ] Because of the similarities of these waves with those 
at the bow shock and the likelihood that this magnetic struc- 
ture originated at the bow shock, we suggest that both sets 
of waves originate from the same free energy source. This is 
discussed in the following section. 


4. Discussion 

[ 23 ] We have presented a survey of STEREO and Wind 
burst captures in association with magnetic field gradients 
at the bow shock and downstream in the magnetosheath. 
Burst captures at the bow shock consist of three types of pre- 
viously identified waves, ion acoustic waves, electrostatic 
solitary waves, and whistler mode waves, as well as newly 
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Figure 4. Comparison of Wind burst captures of waves with cyclotron harmonics (a) in the nightside 
magnetosphere at GSE=[-3,-6,-0.5] (R E ) and (b) at a bow shock crossing. 


identified electron cyclotron harmonic waves. Cyclotron 
harmonic waves were also observed just downstream of the 
bow shock at magnetic field gradients in the magnetosheath. 
These waves are easily identified by their distinctive comma- 
shaped polarizations in the plane containing the magnetic 
field and maximum variance electric field. The variety of 
bow shock crossing locations (and hence shock normal 
angles) with these cyclotron harmonic waveforms (Table 1) 
suggests that they may be ubiquitous at the bow shock and 
also within the magnetosheath (Table 2). 

4.1. Comparison With Traditional Bernstein Waves 

[ 24 ] The waveforms we have presented show signs of 
cyclotron harmonics yet they differ from typical Bernstein 
waves of the sort observed within the magnetosphere in a 
couple of important ways. Figure 4 shows a Wind burst com- 
parison between a more typical Bernstein wave, observed 
in the nightside magnetosphere and a cyclotron harmonic 
wave observed at the bow shock. Some differences are 
immediately apparent. FFTs show that both waves have 
harmonically spaced peaks in power. The traditional Bern- 
stein wave in Figure 4a has power strongly focused at these 
peaks with little power at other frequencies. In this particu- 
lar example a majority of the total wave power is near f ce , 
giving the waveform a sinusoidal quality. The bow shock 
FFTs, in contrast, have significant wave power off-peak as is 
seen in Figure 4b. These waveforms typically have a choppy, 
asymmetrical quality. 

[ 25 ] Another difference is seen by comparing polariza- 
tions. Traditional Bernstein waves tend to show simple 
polarization, primarily perpendicular to the magnetic field. 
The bow shock cyclotron harmonic waveforms, in con- 
trast, often have a complex comma-shaped polarization with 
significant power both parallel and perpendicular to the mag- 
netic field (e.g., Figures Id and 3e). Since harmonics tend 


to be more prominent in the directions perpendicular to 
the magnetic field, this suggests an interaction between two 
waves, a perpendicular Bernstein wave, and a parallel wave 
without harmonics. One likely candidate for the parallel 
wave is ion acoustic waves which have been shown to be 
ubiquitous at the bow shock. 

[ 26 ] These observations suggest a fundamental difference 
between traditional Bernstein waves and the waves observed 
at the bow shock. The differences may be due to the gen- 
eration mechanism. Traditional Bernstein waves are known 
to grow from loss cone distributions [Horne and Thome, 
2000] while the bow shock waves, given their close associa- 
tion with magnetic field gradients, likely get their free energy 
from beam sources. 

[ 27 ] A possible generation mechanism is the electron 
cyclotron drift instability (ECDI) [Forslund et ah, 1970] 
which involves the coupling of Bernstein and ion beam 
modes that occurs when, for example, shock-incident elec- 
trons see (in their rest frame) ion acoustic waves from 
shock-reflected ions Doppler-shifted to near f ce . This can 
result in a strong interaction between electron and ion 
species. A recent particle-in-cell simulation by Muschietti 
and Lembege [2013], designed for conditions that allow the 
growth of small-scale waves (realistic values of the proton to 
electron mass ratio, the plasma to electron cyclotron fre- 
quency ratio, etc.), has indicated that the ECDI may be much 
more important than previously thought at the bow shock 
transition region. The waveforms in this paper are consis- 
tent with those in Muschietti and Lembege [2013]. Namely, 
the estimated wavelength of a few times 100 m is consistent 
with that predicted (A ~ 60 A/,), power is observed to peak in 
close association with f ce , and the comma-shaped polariza- 
tions indicate the presence of both Bernstein waves and ion 
acoustic waves. Thus, the ECDI represents a possible source 
for these waves. 
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[ 28 ] The simulations by Muschietti and Lembege [2013] 
indicated that the ECDI has an inverse cascade of harmon- 
ics with time, i.e., the frequency range of unstable harmonics 
decreases with time. Our observations indicate that the lower 
harmonics tend to dominate over higher harmonics, suggest- 
ing that large amplitude ECDI waves are formed during later 
stages of ECDI evolution. 

[ 29 ] This study represents the first observational evidence 
that the ECDI may occur at the bow shock. The only prior 
observation was presented by Wilson et al. [2010], who 
reported ECDI waves from Wind burst data at an inter- 
planetary shock. Future validation of the ECDI at the bow 
shock and magnetosheath is important because, as pointed 
out by Muschietti and Lembege [2013], ECDI waves, though 
they carry only a small fraction of the energy density 
of the electron and ion beams, are expected to be effec- 
tive at transferring energy and momentum between the 
species, ultimately resulting in energy dissipation. True val- 
idation of the ECDI may be addressed from a more detailed 
analysis of the Wind burst data set in conjunction with 
particle data. 

4.2. Identification of Cyclotron Harmonic Waves 

[ 30 ] We conclude with a discussion of the identification of 
electron cyclotron waves at the bow shock. Our study is sim- 
ilar in many ways to Hull et al. [2006] who analyzed high 
time resolution waveform data from Polar at the bow shock 
during a time of extremely high solar wind velocity. Their 
analysis identified waves from 100 to 4000 Hz with ampli- 
tudes of up to 80 mV/m, lasted from 10 to 30 cycles, and 
were observed at magnetic field gradients. An interferomet- 
ric analysis showed phase velocities consistent with the ion 
acoustic speed and wavelengths of roughly half the electron 
cyclotron radius. The authors thus concluded that the wave- 
forms were likely ion acoustic. This result is consistent with 
past studies that identified ion acoustic waves as the primary 
component of wave power at >200 Hz near the bow shock 
ramp [ Rodriguez and Gurnett, 1975; Wygant et al., 1987] 
and at interplanetary shocks \ Wilson et al., 2007]. 

[ 31 ] We note, however, that waveform properties deter- 
mined by Hull et al. [2006] are also consistent with those 
of cyclotron harmonic waves observed in this study. When 
harmonics are not apparent and polarization is not pre- 
dominantly along the magnetic field, it can be difficult to 
distinguish between cyclotron harmonic waves and Doppler- 
shifted ion acoustic waves, though power spectra above f ce 
tend to be broader for the cyclotron harmonic waves than 
for traditional ion acoustic waves [e.g., Wilson et al., 2010, 
Figure 10] and they typically exhibit the comma-shaped 
polarization discussed previously. More sophisticated analy- 
sis techniques [e.g., Giagkiozis et al., 2011] may be needed 
to identify the difference between these two types of waves 
when clear harmonics are not present. 

[ 32 ] We now discuss reasons why the unambiguous iden- 
tification of harmonics may be difficult even when the burst 
captures in fact contain cyclotron harmonic waves. The 
first difficulty is that, as discussed previously, the instabil- 
ity that generates these waves may involve a combination of 
cyclotron harmonic waves and ion acoustic waves, for exam- 
ple, in the ECDI. Ion acoustic waves are often not polarized 
strictly along the magnetic field [e.g., Fuselier and Gurnett, 
1984; Hull et al., 2006] and may make a significant 


contribution to the perpendicular wave power, tending to 
mask harmonics which may otherwise be distinct. 

[ 33 ] A second difficulty involves details of which harmon- 
ics are unstable. In both the traditional Bernstein instability 
and the ECDI unstable frequencies can, depending on con- 
ditions, fall anywhere from nf ce to (« + 1 )f ce where n is an 
integer, and successive harmonics are not necessarily mul- 
tiples of the first harmonic. As an example, Muschietti and 
Lembege [2013, Figure 2] show that the frequency of each 
unstable harmonic for the ECDI depends on the speed of 
the shock-reflected ions, and the unstable frequencies will 
not be exact harmonics of f ce . This often does not present an 
obstacle toward identifying traditional Bernstein waves (as 
in Figure 4a) because the harmonics can dominate the spec- 
trum but can be for the cyclotron harmonic waves observed 
at the bow shock (as in Figure 4b) because the harmonics do 
not dominate the spectrum to the same extent. 

[ 34 ] Finally, even if the unstable frequencies are exact 
multiples of the fundamental and are confined to the per- 
pendicular plane, we still may not be able to observe 
them for the following reason. Cyclotron harmonic waves 
tend to be short wavelength and thus can be significantly 
Doppler-shifted by the solar wind flow past the observing 
spacecraft. When these waveforms propagate with a compo- 
nent along the solar wind velocity, higher harmonics, which 
have shorter wavelengths, can be Doppler-shifted more than 
lower harmonics. This effect will tend to mask the presence 
of harmonics. If there is to be any hope of observing a pre- 
served harmonic spacing in power, the wave vector must be 
largely perpendicular to the solar wind velocity. This is the 
case with the examples presented in this paper but is not 
expected to be the case in general. The above reasons show 
that we are able to detect the presence of cyclotron harmon- 
ics in waves with comma-shaped polarizations only under 
special circumstances. 

[ 35 ] Waves with comma-shaped polarizations are the 
most common type of wave seen during the STEREO bow 
shock crossings, regardless of bow shock geometry, which 
suggests that they are a common feature of the bow shock 
transition region. Furthermore, ion acoustic waves at the 
bow shock may often be coupled to Bernstein waves via 
the ECDI. 

5. Conclusions 

[ 36 ] We have presented the first observations of electron 
cyclotron harmonic waves at the Earth’s bow shock. The 
waves are large amplitude (up to >500 mV/m peak-peak) 
and are observed at magnetic field gradients at nearly every 
bow shock crossing, regardless of the shock normal angle. 
These waves are also observed in the magnetosheath in 
association with a magnetic field structure, likely convected 
downstream from the bow shock. 

[ 37 ] The cyclotron harmonics are often most prominent 
in the directions perpendicular to the magnetic field. FFTs 
of Wind burst captures with electron cyclotron harmonics 
often show a similar response on both the short (15 m) 
and long (100 m) spin plane antenna at low harmonics, but 
with a decreasing response on the long antenna at higher 
harmonics. This indicates that the first harmonic has a wave- 
length of a few times 100 m, consistent with the electron 
cyclotron radius. 
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[38] The burst captures generally have significant power 
along the magnetic field direction as well as perpendicular 
to it, resulting in a distinctive comma-shaped polarization. 
These observations suggest that the waves are a combination 
of a perpendicular Bernstein wave and a parallel ion acoustic 
wave. The most likely mechanism for the creation of these 
waves is the electron cyclotron drift instability which has 
recently been predicted to occur at the bow shock. 

[ 39 ] Our observations show that a majority of burst wave- 
form captures observed by STEREO during the bow shock 
crossings exhibit this comma-shaped polarization, outnum- 
bering other types of high-frequency waves previously iden- 
tified like whistler mode waves, ion acoustic waves, and 
electrostatic solitary waves. Therefore, we expect that they 
are an important component of the high-frequency electro- 
static spectrum at the bow shock. 


Appendix A: TDS Signal Calibration 

[ 40 ] In this appendix we discuss the calibration of the TDS 
electric field data, first on STEREO and then on Wind. 

[ 41 ] The measurement of external electric field signals 
by the TDS on each STEREO spacecraft involves the cou- 
pling of the antenna to the signal and the spacecraft body, 
and the modification of the signal as it proceeds through the 
physical circuitry of the preamplifier and the TDS instru- 
ment with its 16 bit analog-to-digital converter. Ignoring 
for the moment details of the antenna coupling, a signal 
measured by an antenna arrives at the front end of the pream- 
plifier as the time-dependent voltage V in (t). It then passes 
through the preamplifier and is then converted to counts 
as a function of time N(t) by the analog-digital converter 
in the TDS. This signal is related to the input voltage via 
the equation V m (t) = N(t ) ( ±2 jj/ lllls j, where the last term is 
the resolution of the 16-bit analog-to-digital converter. The 
frequency-dependent gains of the preamplifier and the TDS 
circuitry are then applied to the Fourier transform of the 
signal via FFT( V m (t)) ^ — L— j and the result is inverse 
transformed back to the time domain giving V out . This proce- 
dure is applied to the signal on each antenna. This calibration 
is straightforward and depends only on the characteristics of 
the physical circuitry. 

[ 42 ] The conversion of V out to an electric field on each 
antenna is given by the equation E oat = V out /L eS , where L eS 
is the effective antenna length, also known as the “elec- 
trical length.” The effective length differs from the phys- 
ical antenna length due to shorting effects from the finite 
antenna capacitance C A , antenna coupling to the plasma 
sheath resistance Rs and capacitance C$, antenna coupling 
to the antenna enclosures via the base resistance Rg and 
capacitance Cg, and the capacitance of the antenna cable 
and preamp input (bundled up with Cg in the term C str ay) 
[Bale et al., 2008, Figure 2]. A proper measurement of a 
three-dimensional electric field signal requires knowledge of 
the frequency-dependent effective lengths of each antenna 
to determine wave amplitude, and the frequency-dependent 
effective angles between the antennas which allows the rota- 
tion of the data from the antenna coordinate system into 
the orthogonal spacecraft coordinate system via the equation 
below Figure 14 of Bale et al. [2008]. 


[ 43 ] At high frequencies / » f pe , where f pe is the electron 
plasma frequency, the antennas are in the capacitively cou- 
pled regime where they are weakly coupled to the plasma 
and the antenna gain is described by G ~ C A /(C A + C stray ) ~ 
1/2. The effective lengths and angles at these frequencies 
are constant over a large frequency range and are well 
determined [Bale et al., 2008, Figure 14]. 

[ 44 ] At lower frequencies an external field is coupled to 
the antennas through the plasma sheath resistance Rs and 
capacitance C$. At very low frequencies (less than a few 
hundred Hz) the system enters the resistively coupled regime 
where Cs can be ignored but Rs is important. The gain in 
this regime is given by G ~ Rg/(Rg + Rs). Rs is a sensitive 
function of the electron thermal current but its effect can be 
minimized by maximizing the base resistance. This was not 
done on STEREO, however, because the accurate measure- 
ment of low-frequency electric fields was not a high mission 
priority. The fundamental frequency of the harmonic waves 
presented in this paper ranges from a few hundred Hz to a 
few kHz, near the estimated transition from the capacitively 
to resistively coupled regimes. We have determined that the 
variation in the effective antenna lengths and angles from 
their capacitively coupled values at these frequencies is min- 
imal for the observed waves, and thus, we have used the 
effective antenna lengths and angles defined in Bale et al. 
[2008, Figure 14], 

[ 45 ] Contrary to STEREO, the Wind antennas measure the 
cyclotron harmonic waveforms in the capacitative coupling 
regime where the effective lengths are well determined. 
From the start of the mission until 3 August 2000, the effec- 
tive lengths were 43 m for the long (100 m physical length) 
tip-tip spin plane antenna and 4.5 m for the short (15 m 
physical length) tip-tip spin plane antenna. However, the 
long antenna was partially sheared off by a dust impact on 
3 August 2000, reducing its effective length to 28 m. This 
occurred again on 25 September 2002, further reducing its 
effective length to 25 m. The effective length of the short 
spin plane antenna is unchanged from its original value. The 
Wind TDS observations presented in this paper are from 14 
November 2001, and thus, we used effective lengths of 28 m 
and 4.5 m. 

Appendix B: Instrumental Response 
to Short Wavelengths 

[ 46 ] Since it is possible for harmonics to be an artificial 
instrument response we now show that the harmonics are 
indeed a real phenomenon. 

[ 47 ] There are two processes in which an antenna/receiver 
can artificially produce harmonics in response to an input 
signal. The first involves electronic nonlinearities created 
when an antenna responds to a large amplitude input sig- 
nal. The antennas on STEREO and Wind were designed to 
respond linearly to large amplitude waves of the type pre- 
sented in this paper. It is unlikely that either antenna system 
would create harmonics at the mV/m level in response to 
large amplitude (but unclipped) waveforms. 

[ 48 ] The second process involves a harmonically varying 
antenna response to waves with wavelengths on the order 
of or shorter than the effective antenna length L e g, as dis- 
cussed in Appendix A. The antenna response at the preamp 
input is 8 VJE 0 L eS , w lie re E„ is the applied electric field. The 
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response is a function of the normalized wave number x = 
kL e ff / 2 and is plotted in Gurnett [1998, Figure 9] for a cylin- 
drical dipole. Note that the Wind antennas are cylindrical 
dipoles but the STEREO antennas are cylindrical monopoles 
with the spacecraft body acting as a grounded plane. Bale 
et al. [2008] showed that this is a good approximation to a 
dipole if the antenna is grounded to the spacecraft body by 
a base capacitance Cb of 32 pF. This effect is accounted for 
implicitly via the usage of the effective lengths determined 
by Bale et al. [2008] in calibrating the burst waveforms. Val- 
ues of the normalized wave number less than unity result 
in a continuously decreasing antenna response as a func- 
tion of frequency while values greater than unity result in a 
harmonically varying antenna response. 

[ 49 ] We now calculate the normalized wave number x by 
assuming that the wavelengths of the cyclotron harmonic 
waves are roughly the electron cyclotron radius (~ 500 m), 
consistent with the observed antenna response on Wind. This 
gives us a wave number of k ~ 12 knr 1 . For the short 
STEREO antennas (Z, e ff ~ 1 .2 m) the normalized wave num- 
ber is x < 0.01 « 1 and the antenna response to this wave 
will be a smoothly decreasing function of frequency. For the 
longer Wind spin plane antenna (Z. eft - = 28 m) and shorter 
spin plane antenna (L C |f = 4.5 m), the normalized wave num- 
bers are x ~ 0.2 and x ~ 0.03, respectively. Thus, the 
Wind antennas, like the STEREO antennas, will respond as a 
smoothly decreasing function of frequency. Higher harmon- 
ics have larger wave numbers and if we extend this analysis 
out to the third harmonic, then we have a wave number of 
k = 2? ■ 12 km 1 and an effective wave number of v ~ 0.2 for 
the short antenna pair which will still result in a smoothly 
varying frequency response. For wavelengths much shorter 
than 500 m, the response on long Wind spin plane antenna 
may exceed unity. However, the effective lengths of the 
two Wind spin plane antennas differ by a factor of nearly 
7. This means that the harmonic response of each antenna 
will be different and it is unlikely that both antennas would 
see the same harmonics, as is typically observed, if they 
were artificial. 

[ 50 ] Finally, we have shown that the harmonics are more 
likely to be observed in the directions perpendicular to the 
magnetic field than parallel to it. Artificial harmonics are 
unlikely to be organized in this manner. The above consid- 
erations indicate that the harmonics are innate to the burst 
waveforms and not artificially created. 
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